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Abstract 

We present a lattice QCD calculation of phase shift including the chiral 
and continuum extrapolations in two-flavor QCD. The calculation is carried 
out for 1 = 2 S-wave vrvr scattering. The phase shift is evaluated for two 
momentum systems, the center of mass and laboratory systems, by using 
the finite volume method proposed by Liischer in the center of mass system 
and its extension to general systems by Rummukainen and Gottlieb. The 
measurements are made at three different bare couplings /3 = 1.80, 1.95 and 
2.10 using a renormalization group improved gauge and a tadpole improved 
clover fermion action, and employing a set of configurations generated for 
hadron spectroscopy in our previous work. The illustrative values we obtain 
for the phase shift in the continuum limit are (^(deg.) = —3.50(64), —9.5(30) 
and —16.9(64) for ■^^(GeV) = 0.4, 0.6 and 0.8, which are consistent with 
experiment. 

PACS number(s): 12.38.Gc, 11.15.Ha 
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I. INTRODUCTION 



Calculation of scattering phase shift is an important step for expanding our understand- 
ing of the strong interaction based on lattice QCD beyond the hadron mass spectrum. For 
scattering lengths, which are the threshold values of phase shifts, several studies have already 
been carried out. For the simplest case of the two-pion system, 1 = 2 S-wave scattering 
length has been calculated in detail [1-8] including the continuum extrapolation [6,7]. There 
is also a pioneering attempt at / = scattering length [3], which is much more difficult due 
to the presence of contributions from disconnected diagrams. 

The first pioneering study of / = 2 S-wave scattering phase shift was made by Fiebig 
et al. [9]. They estimated the phase shift from the two-pion effective potential calculated 
in the lattice simulations. A direct calculation of phase shift on the lattice without re- 
course to efi'ective potentials is possible if one uses the finite volume method proposed by 
Liischer [10,11]. In this method phase shift is obtained from the two-pion energy eigenval- 
ues for finite volume. However, it is non-trivial to extract the energy eigenvalues from the 
time behavior of the two-pion correlation functions, because the correlation functions have 
multi-exponential time behaviors due to presence of multiple states with the same quantum 
numbers. Recently we presented a direct calculation of the phase shift [8] in quenched QCD. 
In this calculation the diagonalization method proposed by Liischer and Wolff [12] was used 
to solve the problem of the multi-exponential behavior. Very recently Kim reported on his 
preliminary results on G- and i7-period boundary lattices, where the problem can be avoided 
by the boundary condition [13]. 

In all previous studies of the scattering phase shift calculations were carried out with 
quenched approximation nor the continuum limit was taken. It is predicted in chiral per- 
turbation theory that unphysical chiral divergences appear in the scattering length in the 
quenched theory due to the lack of unitarity [14,15]. The unphysical divergences can also 
occur in the phase shift. While the presence of such pathologies has not been numerically 
confirmed in actual lattice calculations with quenched approximation, we should make our 
study in full QCD in order to avoid such uncontrollable quenching problems. 

In this paper we present a calculation of the physical scattering phase shift for 7 = 2 
S-wave two-pion system including the dynamical quark effects and taking the continuum 
limit. We use the full QCD configurations previously generated for a study of light hadron 
spectrum [16] with a rcnormalization group improved gauge action and clover fcrmion action 
with a tadpole improved clover coefficient. The phase shift is evaluated by the finite volume 
method as in the previous work [8]. In order to obtain the phase shift at several energies 
from one full QCD configuration, the calculations are carried out for two types of momentum 
systems. One of them is the center of mass system, where the total momentum of the two- 
pion system is zero. In the other system where the total momentum is fixed to a non-zero 
value, we use the method proposed by Rummukainen and Gottlieb [17] which is a simple 
extension of that by Liischer in the center of mass system to general momentum systems. 
We shall refer to this system as laboratory system in this paper. 

This paper is organized as follows. In Sec. II we describe the finite volume method 
for both momentum systems for evaluating the phase shift from the two-pion energy. We 
also explain the diagonalization method to extract the two-pion energy eigenvalues. The 
parameters of the calculation in this work are given in Sec. III. In Sec. IV we show data 
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for the pion four-point functions and the effect of diagonahzations for the two momentum 
systems. We then present the results of the scattering length and the phase shift at each 
and in the continuum limit. In Sec. V we briefly summarize this work. 



II. METHODS 
A. Finite volume method 

1. Center of mass system 

In the center of mass system the total momentum of the two-pion system is zero. The 
energy eigenvalue of two pions in a finite periodic box without two-pion interactions is 
given by 



E,, = 2^ml+pl pi = {271/ LYn, neZ, (1) 

In the interacting case the total momentum is also zero and the energy eigenvalue of the 
n-th state is given by 

E„ = 2^m2+p2, p2^(27r/L)%, n ^ Z. (2) 

The energy eigenvalue is written as that of the free two-pion case with momentum and 
— p„, but the quantity n — (L/27r)^ • is no longer an integer. Liischer [10,11] found that 
the momentum satisfies the relation 

tan(5(pj = — — — , (3) 

ZoQ[l]n) 

where is the S-wave scattering phase shift in the infinite volume and 

Zoo(5;n) = 4= E(^'-^)"'- (4) 



„6Z3 

The calculation method of Zoo(s; ft) is discussed in Appendix A. Using Eq.(3), we can obtain 
the phase shift from the energy eigenvalue calculated in the lattice simulations. 



2. Laboratory system 

Let us consider a two-pion system with a non-zero total momentum P 7^ in a periodic 
box L^. We shall refer to this system as the laboratory system in the following. In the 
laboratory system the energy eigenvalue for the n-th energy state without interaction is 
given by 

= v(2v^)p^:„ez^ L/(27r)p^„ e z^ (s) 

where pf^ is the i-th pion momentum of the n-th energy state in the laboratory system, 
which takes discrete values due to the periodic boundary condition in finite volume. The 
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two-pion interaction shifts to E^^ as in the center of mass system. By the Lorentz 
transformation with a boost factor 



7 = </v(<)^ - p^ (6) 

the energy in the center of mass system can be obtained as 

En = Elh. (7) 

It is noted that the total momentum P is not shifted by the two-pion interaction. 
The center of mass momentum is determined from En by 



En^2^ml+pl, pI = {271 /Lfm, m ^ Z, (8) 

where fn is also not an integer. Rummukainen and Gottlieb [17] found that the momentum 
Pn is related to the phase shift in the infinite volume through the relation 



3/2 /= 



m 



where 

^oo(«;^)--i= E(^'-^)"^- (10) 

The summation for r is carried out over the set 

^{r\r^^-\n + d/2),ne Z^} (11) 

where d = L/(27r)P. The operation 7"^ is the inverse Lorentz transformation : 7~^n = 
1/7 ■ n|| + n_L , where n|| = (n ■ d)d/d'^ is the parallel component and = n — n|| the 
perpendicular component of n in the direction d. The calculation method of Z^Q^sifn) is 
discussed in Appendix A. Substituting d = and 7 = 1 into Eq.(9), we obtain the formula 
for the center of mass system Eq.(3). 



B. Extraction of energy eigenvalues of the two-pion system 

In order to extract the two-pion energy eigenvalues in the system with a total momentum 
P, we construct the pion four-point function by 

where we omit the index for the total momentum. The operator Qn{t) is a two-pion operator 
at time t for the n-th energy eigenstate which is defined by 

^n{t)= E vr(pi,t)7r(p2,t)/[ E 1]' (13) 

Pl,P2GS'n Pl,P2SS'n 
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where 7r(p, t) is the pion operator at time t with momentum p. The pi and p2 are the 
momenta on the lattice and take discrete vahies. i.e., L/(27r)pi £ Z^, L/(27r)p2 G Z^. The 
summation for these momenta is taken over the set 



^n = {Pl,P2| V< + P? + V< + Pi = ^r, Pl + P2 = P}, (14) 

where is fixed at the energy of the n-th energy state with the total momentum P in the 
free two-pion case. This summation of momenta projects out the A'^ representation of the 
rotation group on the lattice, which equals the S-wave states in the continuum, ignoring the 
states with higher angular momentum. For example, states with angular momentum I > 4 
are ignored for L/{2tt)P = (0, 0, 0), and / > 2 for L/{2tt)P = (1, 0, 0) and (1, 1, 0). 
For the source we use a difi^erent operator defined by 

1 Nr 

^^f^n^) = T^EMPi,n,t,e,)vr(p2,n,t,r/,), (15) 

.7=1 



where 



<p,t,Q = ^[E^(^'0e^'''^ej(x)l75[E?(y'0e.(y) 



(16) 



Here pi_„, p2,n are fixed to one of the values in the set 5'„ for the sink operator Jln(i) defined 
in Eq.(13). The functions ^j(x) and r)j{x) are complex random numbers in three-dimensional 
space, whose property is 

The pion propagator is also calculated with the random number as 

1 Nr 

Gl^'^'^Kt) = — E(0k(p„,t)7r(-p„,t5,e.)|0). (18) 

When the number of random noise sources is taken large or the number of gauge 
configurations are large, we expect 

G<''-\t) ^ Glit) = (0|7r(p„,t)7r(-p„,t5)|0) (19) 

and the four-point function will be symmetric under the exchange of the sink and the source 
indices n and m. We fix Nf> = 2 in all calculations. The number of the gauge configurations 
is from 380 to 725 depending on the lattice spacing as shown in Sec. III. 

The four-point function can be written in terms of the energy eigenstates with the 
total momentum P as follows: 

Gnmit) = E ynjVjme-'^^^''''^ V^m = (H,- |O^(0) |0) , (20) 
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where is the energy eigenvalues with the two-pion interaction. Since the matrix Vj^n is 
not diagonal generally, the four-point function contains many exponential terms and is not 
diagonal with respect to n and m. 

A simple method for extracting the energy eigenvalues is diagonalization of G{t) at each 

t. Liischer and Wolff [12] found that the eigenvalue is given by A„(t) = exp(— ii^^t) x {1 + 

p p p 

(9(exp[— A„t])} with A„ = min^^ml-E^ — E^\. In order to extract from the eigenvalue 

by a single exponential fit, we have to analyze G{t) in the large t region where 0(exp[— A„i]) 

terms can be neglected. However, it may be very difficult to employ the fitting with this 

method due to loss of statistics of G(t) in the large t region. 

Liischer and Wolff also proposed another method [12], which is a diagonalization of the 

matrix M(t,to) constructed from G{t) by 

M{t, to) = G{to)-'/'G{t)G{to)-'/', (21) 

where to is a reference time. The eigenvalues of M(t, to) equals 

A„(i,io) = e-^"(*-*°). (22) 

without 0(exp[— A„t]) terms. Therefore after the diagonalization of M(t,to) we can extract 
the two-pion energy eigenvalue by a single exponential fitting of Xn{t,to). In this work we 
adopt the second method. 

In the diagonalization method we assume that the two-pion states dominate and effects 
from other states can be neglected in G{t) for the considered range of t. For example, 
our analysis loses its validity when the center of mass two-pion energy is over the inelastic 
scattering limit, e.g., E > 4m^. In this case G{t) is dominated by the four-pion ground 
state rather than the two-pion state in the large t region. Further we should pay attention 
to the states contained with the excited pion state, such as 7r7r*(1300). The existence of such 
undesirable states is discussed in Sec. 111. 

Since we cannot calculate all components of G{t) in the actual lattice calculations, a 
cut-off of the state N has to be introduced. We expect that the components of Gnm{t) 
with n,m < k dominate for the k-ih eigenvalue Xk{t) in the large t and to region, while 
the components n.m > k are less important. We set t and to large and investigate the 
dependence for the cut-off N > k. 



III. PARAMETERS 



We calculate the scattering phase shift on the gauge configurations previously generated 
including two flavors of dynamical quark effects for the study of the light hadron spec- 
trum [16]. This work employed a renormalization group improved gauge action and clover 
fermion action with a tadpole improved clover coefficient Csw, which we also use in the 
present study. The gauge action is constructed in terms of the 1x1 and 1x2 Wilson loops 
W^""^ and W^""^, 

Sg = ^ f CO E K^^'i^) + E K^^'i^)] ■ (23) 
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The coefficient ci is fixed to ci = —0.331 by an approximate renormafization group analysis 
and Co = 1 — 8ci = 3.648 by the normahzation condition [24] . The bare gauge couphng f3 is 
defined by f3 = Q/g^. For clover fermion action [25] a mean-field improved clover coefficient 
csw = (PF^^^)"^/^ is adopted, where the plaquette 1^^^^ = 1 - 0.8412//3 is determined by 
one-loop perturbation theory [24]. 

The parameters for the configuration generation are summarized in Table I. The con- 
figurations were generated at j3 = 1.80 and 1.95 and 2.10, and for four hopping parameters 
K corresponding to mj^/mp a; 0.6 — 0.8 in each (3. The lattice spacing is estimated from 
the p meson mass and equals a — 0.2150(22), 0.1555(17) and 0.1076(13) fm, respectively. 
The lattice size x T at each /3 is 12^ x 24, 16^ x 32 and 24^ x 48, which correspond to 
a 2.5^ fm^ lattice. The periodic boundary conditions are imposed both in the spatial and 
time directions. 

The quark propagators are calculated with the periodic boundary condition in the spatial 
directions, and the Dirichlet boundary condition in the temporal direction. The source 
operator Qj^R^ts) is set at ^5 = 4 for /3 = 1.80 and /3 = 1.95, and ^5 = 6 for /3 = 2.10 to 
reduce effects from the temporal boundary. In order to avoid effects from excited states the 
reference time for the diagonalization introduced in Eq.(21) is fixed to large value; to = 10, 
12 and 16 at /5 = 1.80, 1.95 and 2.10, respectively. At (3 = 1.80 for the lightest quark mass 
we carry out extra measurements to reduce statistical errors in which the source operator 
is located at ts + T/2 and the Dirichlet boundary condition in the temporal direction is 
imposed at T/2. We average over the two measurements for the analysis of the pion four- 
point functions and the pion propagator. 

In order to extract the phase shift at various momenta from a single full QCD configura- 
tion, the calculations are carried out in three momentum systems, the center of mass (CM) 
and two laboratory systems (LI and L2), whose total momenta are 



In Table II we show the momenta chosen for the source operator defined in Eq.(15). 

All elements in Sn, which appear in the summation over momenta in the sink operator 
Q„ defined in Eq.(13), can be obtained from the source momenta by cubic, tetragonal and 
orthorhombic rotations on the lattice for the CM, LI and L2 systems, respectively. 

The center of mass energies of the free two-pion system En in this work are plotted in 
Fig. 1, where the smallest pion mass at /5 = 2.10 is assumed. In all systems the phase 
shift is evaluated at the ground (n = 0) and first excited states (n = 1), which are denoted 
by closed symbols. Other higher energy states (n > 2) plotted by open symbols are used 
to investigate the effects of the cut-off N introduced in the diagonalization. In the two 
laboratory systems, LI and L2, we also calculate the n = 3 state, because the energies of 
these states are very close to the n = 2 state as shown in the figure and effects from these 
states can be comparable. We also plot the location of the inelastic scattering threshold, 
E > 4m7r, and the energy of state 7r7r*(1300), which can appear in the pion four-point 
function and causes loss of validity of the diagonalization method as discussed in Sec. II B. 
We neglect the quark mass dependence of 7r*(1300) in the estimation. As shown in Fig. 1 the 



CM Center of mass system 
LI Laboratory system 
L2 Laboratory system 



L/(27r)P 
(0,0,0) 
(1,0,0) 
(1,1,0) . 



(24) 
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energies of these undesirable states are very far from the measured point of the scattering 
phase shift. 

We estimate the errors of the four-point functions and the pion propagator by the jack- 
knife method. In our study of the hght hadron spectrum [16], we have shown that the 
separation of 50 trajectories covers all the autocorrelations of the configuration, so that in 
the present analysis we also use bins of 50 trajectories in the jackknife method. In the 
actual analysis a bin size of 5 or 10 {j3 = 2.10 at the lightest quark mass) measurements are 
employed, because we skip 10 or 5 trajectories between successive measurements. 



IV. RESULTS 
A. Effect of diagonalization 

In Fig. 2 we plot the absolute value of the pion four-point function Gnm{t) defined 
by Eq.(12) at several values of n and m for the three momentum systems. This figure 
corresponds to m.,^/mp ~ 0.6 at /9 = 2.10. Filled and open symbols indicate positive and 
negative values. In all systems the off-diagonal components of G„to(^) {t^ tu) are not 
negligible compared with the diagonal components (n = m). We also observe that the 
shape of Gnm{t) in the CM and LI systems are very similar. In the L2 system Goo{t) is 
almost the same as Gii(t) for all t. This is attributed to the fact that these energies are 
very close in the free two-pion case. The figure also shows that Gnmii) is symmetric under 
the exchange of the indices as expected. In the following analysis we assume the symmetry 
of Gnm{i) and we average values of the symmetric components. 

In order to investigate the effect of diagonalization we define two ratios Rn{t) and Dn{t) 
as follows, 

where G'^^J^t) is the i-th pion propagator for the n-th energy state, and A„(i,to) is the n-th 
eigenvalue of the matrix M(t, io) defined by Eq.(21) calculated with a finite cut-off of the 
energy state A^. If the pion four-point function contains only a single exponential term, i.e., 
Gnm{t) oc dnm ' exp[— {t — ts)], then the ratio Rn{t) behaves as 

Rn{t) = A ■ exp{-AEl{t - ts)), (27) 

where AE^ = — and A is a constant. If the cut-off of the state N for the diagonal- 
ization is sufficiently large, then the ratio Dn{t) behaves as 

Dnit) = exp(-AE^(t - to)). (28) 

In these cases we can extract the n-th energy shift AE^ from the ratios Rn{t) or Dn{t) by 
a single exponential fitting. 
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First wc focus on the results in the CM system. For the ground state (n = 0) the ratios 
Ro{t) and -Do(^) for all rriT^/mp and (3 are presented in Fig. 3. We divide the ratio -Do(^) 
by Do^ts) for comparison with Ro{t). The cut-off of the state for the diagonalization is 
set at = 2. We also check the cut-off dependence by taking N — 1 and confirm that it 
is neghgible. Prom the figures we find that the diagonalization does not affect the result 
of the ground state and the energy shift can be extracted from the ratio Ro{t) without the 
diagonalization. 

We compare the ratios for the first excited state (n = 1) in the CM system in Fig. 4, 
where the cut-off is set at = 1 and 2. Here we divide -Di(i) by Di{ts) as for n — 0. In 
contrast to the case of the ground state the diagonalization is very effective for the smaller 
quark masses. The ratio Ri{t) for smaller masses rapidly increases, while such behavior 
cannot be seen in Di{t). The cut-ofi' dependence of Di{t) is negligible for whole parameter 
regions as shown in the figure. Hence we can extract the energy shift from the ratio Di{t) 
by a single exponential fitting. 

The ratios in the LI system are plotted in Fig. 5 for the ground state {n — 0); we also 
divide Dn{t) by Dn{ts) as for the CM system. We find that -Ro(^) agrees with Do(t) with 

= 3. We also check the cut-off dependence by taking A^ = 1,2 and confirm that it is 
negligible. This indicates that the energy shift can be extracted without the diagonalization 
as for the ground state in the CM system. For the first excited state (n = 1), however, the 
diagonalization is effective as shown in Fig. 6. We also find that the effect of the cut-off for 
Di(t) is negligible by comparing the results with A^ = 1, 2 and 3. The ratio Di(t) at /3 = 1.95 
and 2.10 for 171,^/171^ ^ 0.6 do not show good exponential behaviors. We consider that this 
is due to insufficient statistics and do not include these data in the following analysis. 

Finally we focus on the L2 system. The behavior of the ratios in this system is essentially 
different from those in the other systems. The ratios Rn{t) and -D„(t) are shown in Fig. 7 
for the ground (n = 0) and in Fig. 8 for the first excited states (n — 1). In these figures 
Dn{t) is divided by Dn{ts)- Since the energies of n = and 1 states are very close in the 
free two-pion case, not only the ratio for n = 1 but also that for n = are affected by the 
diagonaUzation. We also observe that the cut-off dependence is negligible by comparing the 
results with A" = 1, 2 and 3. The ratios Dn{t) for n = 0, 1 for m.,^/mp ^ 0.6 at /3 = 1.95 are 
not clearly exponential in behavior. These results are not included in the following analysis 
considering that these defects are probably caused by insufficient statistics. 

In order to show the effect of diagonalization in the L2 system clearly, we gather both 
ratios Rn{t) and Dn{t) for n = and 1 in Fig. 9. The data at /3 = 2.10 for rriT^/mp ~ 0.6 is 
shown in the figure. Before the diagonalization the two ratios Rn{t) have almost the same 
shape. After the diagonalization, however, the slope of the ratios Dn{t) for n = decreases 
and that for n = 1 increases. This behavior can be understood by considering a simple 
eigenvalue problem for two degenerate states. We assume that the energies for n = and 1 
states in the L2 system have the same energy E in the free two-pion case. Further we neglect 
the effects from higher energy states. This assumption is supported by the independence on 
the cut-off A". In the interacting case the Hamiltonian H for this system can be written as 



where the components of the Hamiltonian are defined by = {^n\H\Vtrn) with the non- 
interacting two-pion state The A and (3 are unknown constants induced by the two- 




(29) 
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pion interaction. The constant A corresponds to the slopes of Rn{t) up to 0{{l3tY). The 
eigenvalues E of H are given by 

E ^ E + A±(3. (30) 

The values of A ± /3 correspond to the slopes of Dn{t) for n = and 1 in the figure. Here 
we can see that the existence of the off-diagonal component of the Hamiltonian /3 causes a 
separation of Dn{t). 

Up to now we have shown that the ratios Dn{t) for any n in all the momentum systems 

— p 

behave as single exponential functions in t and it is possible to extract the energy shift AE'^ 

by a single exponential fitting of the ratios. Our choice of the fitting range and the results 
— p 

of AE^ are summarized in Appendix B. The procedure of calculating the scattering phase 
p 

shift d{p^) from AE^ is the following. First we construct the two-pion energy eigenvalue 
£■„ by E^ = E^ + AE^ , where E^ is the two-pion energy in the free two-pion case. Then 

we evaluate the Lorentz boost factor 7, the center of mass energy E^ and momentum 
P 

from E^ by Eqs.(6), (7) and (8). Finally we obtain the phase shift 5(p,„) by substituting 
j)\ into the finite volume formulae given by Eqs.(3) and (9). The results arc tabulated in 
Appendix B. In the appendix we also quote the 'scattering amplitude' defined by 

tan5(p) E 

A{m^,p)= - ' -iz, (31) 
p z 

where the amplitude is normalized as 

lim A(m7r,p) = ao^T'TT) (32) 
p— >o 

with the scattering length cq. 



B. Result for scattering length 

In the CM system the momentum for the ground state (n = 0) is very small as 
shown in Appendix B. Thus the scattering length oq can be evaluated from the scattering 
amplitude defined in Eq.(31) by A{mT^,p) aoUiT^. In Fig. 10 we plot oo/m^ as a function 
of at each f3, which arc also tabulated in Appendix B denoted by A(m.,r,p)/ml.. A 
significant curvature in the dependence at large (3 is seen. This has not been clearly 
observed in the previous studies of the scattering length. The existence of a large curvature 
renders the chiral extrapolation of ao/m.j^ very difficult. In this work we attempt to fit the 
data with various fitting assumptions to obtain the scattering length at the physical pion 
mass ttItt = 0.14 GeV. 

The prediction of chiral perturbation theory (ChPT) for the dependence has been 
worked out by Gasser and Leutwyler [26]. In one loop order it is given by 



mt 



167rF2 



167r2F2 



L(/.)-Ci-log^ 



1^^ 



(33) 



where F is the pseudoscalar decay constant in the chiral limit, L^fi) is a low energy constant 
at a scale and Cl — T/'2. It is clear that this one- loop formula cannot be naively applied 
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to our results of the scattering length. As shown in Fig. 10 the dependence of ao/rriT, on the 
lattice spacing is comparable to that of m^. This implies that we should extrapolate our 
data with the formula including the 0(a) effect. Such a formula for the pseudo scalar mass 
and the decay constant have been obtained in one-loop order of chiral perturbation theory 
for the Wilson type fermions [29], but that for the scattering length is not yet available. 
Further, even if the 0{a) effect is negligible, it is not clear whether the formula of ChPT is 
applicable for our data calculated in a heavy pion mass region, m,r = 0.5 — 1.1 GcV. Here, 
as a trial, we fit our data with the fitting assumption given by Eq.(33) with three unknown 
parameters F, L{ix) and C/,, and consider the dependence of fitted parameter values on 
the lattice spacing. The results are summarized in Table III, where the scale ji is fixed at 
^ = 1 GeV in the analysis. We find a significant difference between Cl obtained by the 
fitting and that of the prediction from ChPT (Cl = 7/2) at all j3. 

It is expected that the chiral breaking effect of the Wilson type fermion causes a diver- 
gence in the chiral limit as Oo/mTr oc 1/m^ [27]. The large curvature of the scattering length 
may originate from this effect. In Table IV we tabulate the results of the fitting with 

^ = ^ + A,o + A,om^ (34) 

The coefficient of the divergent term Aqq, which comes from the chiral breaking effect, is 
expected to vanish in the continuum limit. However, our results for Aqq increases toward the 
continuum limit as opposed to the expectation. We consider that the effect of chiral breaking 
is not separated from the regular mass dependence. It is a very important future work to 
detect the divergence term through simulations with much higher statistics and closer to the 
chiral limit. We assume that the effect of chiral breaking is small in our simulation points 
in the following analysis. 

Next we attempt to fit our data assuming the following polynomial function in m^, 

= ylio + A2oml + Asoml. (35) 

The results of the fitting are plotted in Fig. 10 and summarized in Table V. At /3 = 2.10 
the value of x^/d.o.f. is large. This indicates that the fitting including higher order terms, 
such as term or higher, is necessary to obtain a more precise value at the physical pion 
mass. Such fitting cannot be carried out in this work, because the number of our simulation 
points is only four. 

In Refs. [2-4] the authors found that the mass dependence of the ratio (/^"*)^ • Oq/^^tt 
was very small, where /^"* is the decay constant measured on the lattice at each m^. In 
Fig. 11 we plot the normalized scattering length defined by 

ao_//n\i?l, (36) 



where /^"* is the decay constants measured in Ref. [16], which are tabulated in Table VI, 
and = 93 MeV. The values of aQ/m^ arc also tabulated in Appendix B; they are written 
under the column for 74(m7r,p)/m^ defined by 



TV 



' £la 

A{m,,p) = [^-^] ■A(m,,p). (37) 
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In the calculation of ao/'"^7r the statistical errors of /|"* are not included, since they are 
small compared with those of the scattering length. We observe in Fig. 11 that 0,^/711^^ is 
almost independent of within the statistical errors. This fact implies a strong correlation 
between the scattering length and the decay constant. The results of a constant fitting for 
ao/m^ are also plotted in the figure and tabulated in Table VII. The values of x^/d.o.f. are 
reasonably small. Especially at = 2.10 it is much smaller than that of the polynomial 
fitting for ao/m^. 

In Fig. 12 we present aom^ at the physical pion mass obtained by the polynomial fitting 
defined by Eq.(35) and ciom^ by the constant fitting as a function of the lattice spacing. We 
also plot the results of the continuum extrapolations, which are summarized in Table VIII. 
The prediction of ChPT [28]: aoTn,r = —0.0444(10) is denoted by the star symbol in the 
figure. We see large 0(a) effects in both flo'^Tr and 0,^11x1^^. Furthermore the 0{a) effect for 
OomTT is opposite to that of aom^ due to large 0(a) effects in /^"*. The difference between 
them decreases going toward the continuum limit. 

As shown in Fig. 12 OomTr in the continuum limit is closer to the prediction of ChPT 
than that of aom„. We should note, however, that the decay constant is introduced in ao/m^ 
to compensate the mass dependence of the scattering length. Extrapolations only with the 
scattering length are difficult as discussed before. 

The decay constant measured in the previous work of Ref. [16], is tabulated in Table VI. 
As shown in the table, the value in the continuum limit seems not consistent with the 
experiment /tt = 93 MeV. One of possible reasons for this discrepancy is an uncertainty 
of the renormalization factor for the axial vector current; in Ref. [16] the perturbative 
renormalization factor was used. This causes a large uncertainty in the prediction of the 
scattering length in the continuum limit. In order to obtain more reliable predictions for the 
physical scattering length, a non-perturbative determination of the renormalization factor 
and higher statics calculations of the scattering length closer to the physical pion mass 
and the continuum limit are needed. Here we present the two values estimated by Oo^tt 
and aonij, = (/^"*//^)^ ■ aomT^ as our results for the scattering length in this work with the 
provisions noted above: 

aom^ = -0.0558(56), (38) 
aom^ = -0.0413(29). (39) 

C. Result for scattering phase shift 

The results for the scattering amplitudes A{m.,:,p) defined by Eq.(31) are plotted in 
Fig. 13. The CM„ in the top left figure refers to the amplitude obtained from the n-th 
energy eigenstate state in the CM system, and the Ll„ and L2„ to those in each laboratory 
system. The ordering of the amplitude, CMq, LIq, L2q, ■ ■ •, Lli, in the other figures is the 
same. The open symbols are excluded from the following analysis. They correspond to the 
data in which we cannot clearly observe a single exponential time behavior of Dn{t) as noted 
in Sec. IV A. The values of A(m^,p) are tabulated in Appendix B. 

We obtain the scattering phase shift 5(p) at the physical pion mass for various momenta 
p in the continuum limit by the following procedure. First we fit the results of 74(m^,p) 
with the following fitting assumption at each lattice spacing. 
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A(m7r,p) = Aioml + ^20"^^ + ^30 + ^oiP^ + ^nm^p^ + A2im^p^. 



(40) 



Then we evaluate the amphtudes at the physical pion mass for various momenta at each 
lattice spacing from the constants A^j obtained by the fitting. Finally the continuum ex- 
trapolation is taken for these amplitudes at each momentum. The results of the fitting with 
the assumption Eq.(40) are plotted in Fig. 13, and the constants Aij are tabulated together 
with x^/d.o.f. in Table IX. The values of x^/d.o.f. are reasonably small at all [3 in contrast 
to the polynomial fit for the scattering length where x^/d.o.f. ai (3 — 2.10 is large as shown 
in Sec. IV B. 

It is possible to make a global polynomial fit with the assumption Eq.(40) for all momen- 
tum systems, so that we do not need to introduce the decay constant measured on the lattice 
as for the scattering length in the previous section. For consistency of the analysis, we also 
analyze A{mn,p) defined by Eq.(37). Here we assume the same fitting assumption Eq.(40), 
but set ^420 = ^30 = A21 — 0, since the mass dependence becomes mild by a compensation 
with that of the decay constant also in this case. The fit results are tabulated in Table X. 

In Fig. 14 both amplitudes at the physical pion mass for various momenta are plotted 
as a function of the lattice spacing, together with the continuum extrapolations. Here we 
choose three momenta, = 0, 0.06 and 0.26 GeV^ which are roughly equal to the momenta 
of CMq, LIq and CMi, respectively. The result at = gives the scattering length aom^. 
We also plot the scattering lengths obtained in the previous section, which arc calculated 
from the data only in CMq. As shown in the figure they are consistent within the statistical 
error at each lattice spacing. There are slight differences, however, in the continuum limit 
between the results obtained from all momentum system at = given by 



and the results from only the CMq data given by Eqs.(38) and (39). We consider that these 
discrepancies represent the uncertainty of the continuum extrapolations using data at only 
three lattice spacings, e.g., the number of degrees of freedom for the extrapolation is one, 
far from the continuum limit. The difference of the two amplitudes tends to vanish in the 
continuum limit. As mentioned in the previous section, however, there is an uncertainty 
in the determination of the decay constant on the lattice. Thus we do not use A{mT^,p) to 
obtain the final results of the phase shift. 

In Fig. 15 the result of the phase shift obtained from the scattering amplitude A{mT^,p) 
in the continuum limit is presented by a dashed line, and associated by a band of error 
bars. The values of the phase shift at several momenta are tabulated in Table XI. Our 
results arc compared with the solid curve [28] estimated with the experimental input, and 
the experimental results [30,31]. The result in the continuum limit agrees with experiment, 
although the errors of our result are large. In order to obtain more precise results simulations 
are needed closer to the chiral and the continuum limits with much higher statistics. 



In this article we have presented our results for the 1=2 S-wave tttt scattering phase shift 
in the continuum limit calculated with two-flavor dynamical quark effects. While errors are 



0.0484(49), 
0.0404(24), 



(41) 
(42) 



V. CONCLUSIONS 
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not small, it is very encouraging to find that the phase shift in the continuum limit shows a 
reasonable agreement with experiment. 

The large errors of our final results arise from the chiral extrapolation and the continuum 
extrapolation. In order to obtain more precise results, simulations are needed closer to the 
chiral and the continuum limits with much higher statistics. Technically, investigating the 
correlation wc found between the scattering length and the decay constant measured on 
lattice, and detecting effects of chiral symmetry breaking with Wilson fermion action are 
important open issues for future work. 

In addition to being a first calculation of realistic scattering quantity based on first 
principles of QCD, the importance of the present work resides in actually showing that 
various technical methods, such as diagonalization of pion four-point functions and use of 
laboratory systems, necessary for practical success of the finite-volume methods work. Thus 
we can envisage a perspective toward extension of the present work for calculations of other 
important scattering process, such as / = 1 and 7 = two-pion systems, systems including 
unstable particles, and scattering with baryons, which are richer in physics content. These 
processes are more difficult from the point of calculations, however, and algorithmic advances 
are probably needed to evaluate complicated diagrams. 
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APPENDIX A: CALCULATION METHOD OF ZETA FUNCTION 

In this appendix we introduce methods for the numerical evaluation of the zeta function 
ZQQ{s;m) defined in Eq.(lO). This problem has already been discussed by Liischer in Ap- 
pendix A of Ref. [10] and in Appendix C of Ref. [11] for the center of mass system, and by 
Rummukainen and Gottlieb in Sec. 5.2 of Ref. [17] for general systems. Our basic idea is 
the same, but our final expression for the zeta function is simpler and more efficient for nu- 
merical evaluations. (We find that there are some typographical errors in the expression of 
the zeta function in Ref. [17]). Very recently Li and Liu have reported a similar calculation 
method of the zeta function in asymmetric box [32] . 

The definition of the zeta function ZQQ(s;m) is 

■ Z^o{s; m)= (r^ - m)-'. (Al) 

The summation for r is carried out over the set 

pd = {r|r = 7-^(n + d/2),ne Z^}. (A2) 

The operation 7"^ is the inverse Lorentz transformation : 7~^n = I/7 • n|| -|- where 
n|| = (n • d)d/ci^ is the parallel component and nj_ = n — n|| the perpendicular component 
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of n in the direction d. The zeta function Z^Q(s;rn) takes a finite vahie for Res > 3/2, 
and ZQQ(l;m), which is used to obtain the scattering phase shift in Eq.(9), is defined by the 
analytic continuation from the region Res > 3/2. 

First we divide the summation in Zqq{s; m) into two parts as 

^(^2 _ -ys ^ ^ (^2 _ -ys ^ ^ (^2 _ -ys ^ ^^3) 
r r'^<m r^>rn 

where the summation over r is carried out with r e P'*. The second term can be written in 
an integral form as follows, 



/' dt f-^e"^ y e-''" - y (r^ - m)-^ + y . (A4) 



r(s) 



-(r^—m) 

T{s)Jo ^" ^ '^{r'^-my 

The second term cancels out the first term in Eq.(A3). Next we rewrite the first term in 
Eq.(A4) by the Poisson's summation formula 

E /(n) = E / d^^/(x)e^^'^"-, (A5) 
and integrating over x yields, 

dt i^-^e^^Ve-*'-' = -2- /' di i^-^e*^ ( -V^' V (_i)n-de-^{7n)Vt. (A6) 
Jo r r s Jo V w ~-:iq 



r(s) Jo 



ez3 



The divergence at s = 1 comes from the n = part of the integrand on the right-hand 
side. We divide the integrand into a divergent part (n = 0) and a finite part (n 7^ 0). The 
divergent part can be evaluated for Res > 3/2 as 

dt i-ie*- ( - ) = V -4 (A7) 

\tj ^s + /-3/2 /! ^ ^ 

The right had side of this equation takes a finite value at s = 1. 

Finally by gathering all terms we obtain the following expression for the zeta function 
at s = 1, 

^ — m 

/•I /7r\3/2 ' 00°^ -77^/2 rrr' 



where Zlnezs is the summation without n = 0. 

Substituting d = and 7 = 1 into the above expression, we obtain the representation of 
the zeta function in the center of mass system appeared in Eq.(3) 

„-(n2-m) „i /7r\3/2 / „ „ ^ ~3/2 
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APPENDIX B: TABLE FOR RESULTS OF SCATTERING LENGTH AND 
SCATTERING PHASE SHIFT IN EACH SYSTEM 



In Tables XII - XVII we tabulate fitting ranges, energy shift AE'^ , center of mass momen- 
tum p^, Lorentz boost factor 7, scattering phase shift S{p), scattering amplitude A{'m.j^,p) 
defined by Eq.(31), and normalized scattering amplitiidc A(mT^,p) defined by Eq.(37) in 
each system for the ground n — and first excited n = 1 states. 
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TABLES 



p — i.oU 


K 


n 1 AC\A 
U. 1404 


r\ ^ A AK 


n 1 /I Qn 
U.i4oU 


n 1 /ino 
U.14Uy 




ml [GeV2] 


W.04 < j 

0.238(1) 


0.571(1) 


U. 1 ooyl J 

0.814(1) 


n ^c\7( 1 ^ 

U.oU I yi- ) 

1.128(1) 


Fitting Range 




10-20 


12 - 20 


12 - 20 


12 - 20 




[xiu LreVJ 


4.Di(^yo j 




n Q A f Q'7\ 

z.o4(^o / j 


O ^ '7f'\ A\ 
Z.l ([14: ) 




[XiU IjCV J 




on Ci( QO^ 
Zy.yyoz, ) 


Z( .\)[o ( ) 


IK 1 / 1 

Zd.i[i ( ) 




[degrees] 


-1.13(34) 


-1.49(22) 


-1.36(25) 


-1.17(11) 






U. -Lc/Ul OO j 


— n 3fi9('35'l 


u.^uui ^y j 


U.^O^l ZU J 


A(m-,r. v„ ) Itt? 

^ J- V " "TT ) Fn / / " "^TT 


[l/GeV^l 

L / ' J 

fl/GeV^l 


-0.82(16) 
-2.47(48) 


1 (^f{(^fi\ 

— i.DO^^lD J 

-0.633(62) 
-2.91(28) 


-0.499(60) 
-3.02(33) 


— o. 1 / ^^-Ly ) 
-0.384(23) 
-2.83(17) 


p — i.yo 




n 1 /1 1 n 

U. i4iU 


n 1 /inn 

U.i4UU 


n 1 Qon 
u.ioyu 


u.io ( 


= 1.268(13) [GeV] 


m2 [GeV2] 


0.582(3) 
0.291(2) 


0.690(1) 
0.573(1) 


0.752(1) 
0.857(1) 


0.804(1) 
1.287(1) 


Fitting Range 




12 - 23 


13-25 


13-25 


13-25 




Twin— 3 r^^AT"! 

[XiU LrCVJ 


i\}.\)0\ (6) 


D.oy^oy ) 


0. {^[oU ) 


O.OD^ZO j 


Pn 


[XiU LreV J 


Afi T/QI \ 

40. / [oi ) 


4i.Z(^4i j 


4i.y(^zz j 


o4.0(zi j 






— Z.OUl^Zo j 


9 1 0/^90^ 


9 1 ^\ 
— Z. lOt^iO j 


— 1. 00^14 j 


^K'lH^Pn) 

/I ( TT),^ T) ^ / 777,^ 


[l/GeV2] 
fl/GeV^l 


-1.195(70) 
-2 63fl5) 


l.O 1 I -LZ J 

-0.759(67) 
-2 39(21) 


Z.Z / I -LU J 

-0.630(29) 
-2 65fl2) 


-0.433(23) 
-2 40fl3) 






U. iooZ 


U.io (4 


U.ioD ( 


U.ioO ( 






u.o loyo) 


n fiQi ( 


U. ( 00(^/ j 


U.oUDl^i ) 




ml [GeV2] 


0.291(1) 


0.605(2) 


0.896(1) 


1.332(2) 


Fitting Range 




18-35 


18-35 


18-35 


18-35 


Tl 

pI 


[xlO-^ GeV] 
[xlO-^ GeV2] 


17.33(54) 

51.3(16) 


10.48(42) 

44.6(17) 


8.49(32) 

44.0(16) 


6.16(30) 

38.9(19) 


^(Pn) 


[degrees] 


-3.17(13) 


-2.62(14) 


-2.57(13) 


-2.17(14) 


Mnn^iPn) 




-0.421(12) 


-0.536(18) 


-0.643(20) 


-7.04(30) 


M^7r,Pn) 

MrnT,,Pn)/ml 

MlT^7r,Pn)/^l 


[l/GeV2] 
[l/GeV2] 


-0.718(20) 
-1.444(38) 
-2.462(66) 


-1.471(49) 
-0.885(31) 
-2.429(85) 


-2.251(72) 
-0.718(23) 
-2.512(82) 


-3.12(13) 
-0.528(23) 
-2.34(10) 



TABLE XII. Results for n = state in the center of mass system CM with energy state 
cut-off N = 2. Two scattering amplitudes are defined by A(m7r,^„) = taii5(p^)/p^ ■ En/2 and 
A{m.jr,Pn) = (/w^V/tt)^ ■ A(m7r,^„), where /^"* is the pseudoscalar decay constant measured on 
lattice and f^, = 93 MeV. 
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/3 _ 1 on 
fj — L.Okj 


K/ 


n I4fi4 


n 144^ 




n 1409 




in^l nip 


f) ^A7(A\ 




u. i ooyi. J 


n RC\7( 1 






f) 9^8fl "1 

L/.^OO I J- ( 


'^7^ (A \ 

XJ.O I -L I -L / 


814('Tl 


1 198fTl 






10 - 18 


12 - 18 


12 - 20 


12 - 20 




1 A lU V J 




1 ^ 9n 9^1 






Pn 


rvin-2 OpV^I 

\/\ J-\J \J\j V J 


24 8n('9fi'l 


94 40('1 


94 380^78 


24 '22'2(T\'\ 






— 14.1(22) 


— 10.7(10) 


—in ^Q((\f,) 


— 9 25(43") 


A \ Tn T) 1 




— n S5'^l'57'> 


— 348l'33') 


—0 38qi'24'> 


—0 387(1 8") 






-1.06(17) 


-1.59(15) 


-2.34(14) 


-2.83(13) 


/3 _ 1 QIX 


K 




IzLOO 


\J.A.Ou\J 


yj.A-O f o 


Cfc — ±.ZrU0l±01 1 Vjrt? V 1 


qrtn 1 nm 
i''"kI ii''p 




U.Uc/UI ± J 










n 9Q1 ('9'> 


U.O 1 Ol X 1 










12 - 23 


13 - 25 


13 - 25 


13 - 25 




[ X lU vjr6 V J 

[xlO-2 GeV^] 




9Q Q/'9n'\ 


99 7r\(Q^\ 






27.55(30) 


27.02(14) 


26.761(79) 


26.389(48) 


"\Pn) 




—20 ^(24) 


—16.7(11) 


— 14 63(63) 


— 11 70(38) 


-'^\'""7TyFn/ 




—0 549 f 68") 


-0 53^38") 


-0 535 (23") 


-0 502(16) 






-1.21(15) 


-1.67(12) 


-2.256(99) 


-2.787(91) 


— 9 1 n 




n 1 '^S9 


n 1 "^74 




n 1 "^^7 


Cfc — ±.0001^^1 1 VjfC V J 


in-j^l nip 

777 2 [npv^i 


yj.o t ui o J 

9Q1 (^ \ 


n fiQi f'^') 

u.uy X I o 1 

fin^i'9') 


n 7'^'^(9') 
8Qfifl "1 

\J.OU\J I J- ( 


n 8nfi( 1 ) 

1 339(9) 


Fitting Range 




18-35 


18-35 


18-35 


18-35 




[xlO-=^ GeV] 


58.8(69) 


44.1(25) 


34.6(14) 


26.80(78) 


pI 


[xlO-2 GeV2] 


25.30(28) 


25.17(13) 


24.969(84) 


24.789(53) 


HPn) 


[degrees] 


-19.8(24) 


-18.7(11) 


-16.99(73) 


-15.44(46) 






-0.530(69) 


-0.626(39) 


-0.654(28) 


-0.697(21) 


Ml^-K,Pn) 




-0.90(11) 


-1.71(10) 


-2.29(10) 


-3.097(94) 



TABLE XIII. Results for n = 1 state in the center of mass system CM with energy state 
cut-off N = 2. Two scattering amplitudes are defined by ^(m^,p„) = tan6{p^)/pj^ ■ En/2 and 
A{m-^,p^) = (/^"*//jr)^ • A{mT^,Pn), where /^"* is the pseudoscalar decay constant measured on 
lattice and f^, = 93 MeV. 
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P = 1.80 


K 


0.1464 


0.1445 


0.1430 


0.1409 


= 0.9176(93) [GeV] 




0.547(4) 


0.694(2) 


0.753(1) 


0.807(1) 




ml [GeV2] 


0.238(1) 


0.571(1) 


0.814(1) 


1.128(1) 


Fitting Range 




10 - 18 


12 - 20 


12 - 20 


12 - 20 




[xlO-3 GeV] 


9.5(10) 


6.72(50) 


4.65(43) 


3.70(19) 


pI 


[xlO-3 GeV^] 


54.16(69) 


58.93(46) 


59.04(45) 


59.56(23) 


7 




1.09436(62) 


1.04481(11) 


1.032523(71) 


1.024028(25) 




[degrees] 


-7.03(76) 


-7.14(52) 


-5.81(52) 


-5.36(27) 






-0.286(30) 


-0.409(28) 


-0.391(34) 


-0.419(20) 


^("^7^,P„) 




-0.862(91) 


-1.88(13) 


-2.35(20) 


-3.07(15) 


P = 1.95 


AC 


0.1410 


0.1400 


0.1390 


0.1375 


= 1.268(13) [GeV] 


■m-^/nip 


0.582(3) 


0.690(1) 


0.752(1) 


0.804(1) 




ml [GeV2] 


0.291(2) 


0.573(1) 


0.857(1) 


1.287(1) 


Fitting Range 




12 - 23 


13 ~ 25 


13 - 25 


13 - 25 




[xlO-=* GeV] 


17.3(11) 


11.99(55) 


9.62(35) 


6.92(30) 




[xlO-3 GeV^] 


61.40(58) 


64.47(36) 


65.73(28) 


65.91(28) 


7 




1.08447(51) 


1.04758(13) 


1.033120(56) 


1.022720(32) 




[degrees] 


-9.29(59) 


-8.55(38) 


-8.20(29) 


-7.12(30) 






-0.392(24) 


-0.473(20) 


-0.540(18) 


-0.566(23) 


HmT:,Pn) 




-0.865(53) 


-1.491(63) 


-2.276(79) 


-3.14(12) 


/3 = 2.10 


K 


0.1382 


0.1374 


0.1367 


0.1357 


= 1.833(22) [GeV] 


m^^/mp 


0.576(3) 


0.691(3) 


0.755(2) 


0.806(1) 




ml [GeV2] 


0.291(1) 


0.605(2) 


0.896(1) 


1.332(2) 


Fitting Range 




18 - 35 


18 - 35 


18 - 35 


18 - 35 




[xlO-3 GeV] 


27.8(12) 


17.7(65) 


14.75(48) 


10.45(38) 


pI 


[xlO-3 GeV2] 


58.78(42) 


61.02(29) 


62.22(26) 


61.97(25) 


7 




1.07874(41) 


1.04214(13) 


1.029504(58) 


1.020373(39) 


KPn) 


[degrees] 


-11.09(46) 


-9.66(33) 


-9.55(29) 


-8.16(28) 


A{m^,Pn) 




-0.478(19) 


-0.563(18) 


-0.660(19) 


-0.680(23) 


A{'m-jT,Pn) 




-0.816(32) 


-1.544(50) 


-2.311(68) 


-3.02(10) 



TABLE XIV. Results for n = state in the laboratory system LI with energy state cut-off 
N = 3. Two scattering amplitudes are defined by A{mTr,Pn) = ^^^S{pJ^)/pJ^ • En/2 and 
A{mTr,Pn) = (/i"*//7r)^ ■ A.{mT^,Pn), where /^"* is the pseudoscalar decay constant measured on 
lattice and f^, = 93 MeV. 
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P = 1.80 




0.1464 


0.1445 


0.1430 


0.1409 


= 0.9176(93) [GeV] 




0.547(4) 


0.694(2) 


0.753(1) 


0.807(1) 




ml [GeV2] 


0.238(1) 


0.571(1) 


0.814(1) 


1.128(1) 


Fitting Range 




10 - 18 


12 - 18 


12 - 20 


12 - 20 




[xlO-3 GeV] 


28.0(63) 


18.1(26) 


14.6(12) 


10.62(59) 


pI 


[xlO-2 GeV^] 


30.65(53) 


30.41(27) 


30.31(14) 


30.061(79) 


7 




1.05169(54) 


1.03245(10) 


1.025510(52) 


1.020014(19) 


HPn) 


[degrees] 


-16.9(39) 


-13.4(20) 


-12.1(10) 


-9.86(56) 


A{m^,Pn) 




-0.405(98) 


-0.406(61) 


-0.412(35) 


-0.379(21) 


^("^7^,P„) 




—1.22(29) 


-1.86(28) 


—2.47(21) 


-2.77(15) 


^ = 1.95 


AC 


0.1410 


0.1400 


0.1390 


0.1375 


= 1.268(13) [GeV] 




0.582(3) 


0.690(1) 


0.752(1) 


0.804(1) 




ml [GeV2] 


0.291(2) 


0.573(1) 


0.857(1) 


1.287(1) 


Fitting Range 




12 - 16 


13 - 20 


13 - 25 


13 - 25 




[xlO-=* GeV] 


57(10) 


34.0(37) 


28.3(16) 


18.98(79) 


pI 


[xlO-2 GeV^] 


34.21(67) 


33.30(28) 


33.26(14) 


32.788(81) 
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1.04789(53) 


1.03371(12) 


1.025787(45) 


1.019069(26) 




[degrees] 


-25.0(47) 


-17.3(19) 


-16.35(99) 


-12.59(53) 


^(m^,p„) 




-0.63(13) 


-0.516(60) 


-0.555(34) 


-0.495(21) 


^("i7r,^„) 




-1.40(29) 


-1.62(19) 


-2.33(14) 


-2.75(11) 


/3 = 2.10 


K 


0.1382 


0.1374 


0.1367 


0.1357 


= 1.833(22) [GeV] 


m^^/mp 


0.576(3) 


0.691(3) 


0.755(2) 


0.806(1) 




ml [GeV2] 


0.291(1) 


0.605(2) 


0.896(1) 


1.332(2) 


Fitting Range 




18 22 


18 - 35 


18 - 35 


18 - 35 




[xlO-3 GeV] 


68(13) 


53.9(48) 


39.3(24) 


32.4(12) 


pI 


[xlO-2 GeV2] 


31.16(58) 


31.23(24) 


30.82(14) 


30.798(89) 
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1.04646(43) 


1.03077(10) 


1.023545(48) 


1.017343(30) 




[degrees] 


-21.8(43) 


-20.7(19) 


-17.1(10) 


-16.32(64) 


^(m^,p„) 




-0.55(12) 


-0.651(63) 


-0.608(39) 


-0.675(27) 


^(?"7r,Pn) 




-0.95(20) 


-1.78(17) 


-2.13(13) 


-3.00(12) 



TABLE XV. Results for n = 1 state in the laboratory system LI with energy state cut-off 
N = 3. Two scattering amplitudes are defined by A{mT^,p^) = tan(5(^„)/^„ • En/2 and 
A{mT^,Pn) = (/i"*//7r)^ ■ ^("^TTjPn)) where /^"* is the pseudoscalar decay constant measured on 
lattice and ^ = 93 MeV. 
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U. -LOU I 


1 ^57 

U. lOO t 


— ^ R'^'\(9'?^ mpvi 
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n fiQi ^'^'1 

\j.\jiy± I o / 


755f9'l 


n 80fiM ) 




777 2 rripV^l 

[v^ev J 


n 9Q1 (^ ) 


n fin5('9'i 


n 8Qfi('i 


1 '^'^9('9) 

± .0(J.^ \ J 


Va tH" t n (T I-< o n rr*^ 




io oU 


lo OO 


io oO 


lo oO 


ft 


GeV] 


17.5(35) 


9.3(15) 


7.83(88) 


4.45(61) 




[xlO-2 GeV2] 


9.48(13) 


10.332(76) 


10.746(50) 


10.932(42) 
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1.13877(72) 


1.07789(27) 


1.05563(10) 


1.039047(74) 


S{Pn) 


[degrees] 


-14.5(41) 


-12.0(33) 


-17.3(56) 


-10.9(40) 


A{m^,Pn) 




-0.52(15) 


-0.55(16) 


-0.95(33) 


-0.70(26) 


M^-K^Pn) 




-0.89(26) 


-1.53(44) 


-3.3(11) 


-3.1(11) 



TABLE XVI. Results for n = state in the laboratory system L2 with energy state cut-off 
-/V = 3. Two scattering amplitudes are defined by A(m7r,p„) = tan(5(^„)/^„ • En/2 and 
A{mT^,Pn) = (/i"*//7r)^ ■ ^("^TTjPn)) where /^"* is the pseudoscalar decay constant measured on 
lattice and ^ = 93 MeV. 
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ft — ^n 


K 


u. ±ooz 


n 1 ^74 


n 1 ^fi7 


n 1 ^f\7 

U. -LOO t 


= 1.833(22) [GeV] 


ml [GeV2] 


0.576(3) 
0.291(1) 


0.691(3) 

0.605(2) 


0.755(2) 

0.896(1) 


0.806(1) 
1.332(2) 


r lining rLctnge 




18 "in 


IS '^1^ 


18 "^1^ 


18 '^1^ 
io 00 




[xin^-^ GpVI 


SI 5('S2'l 


21 2('12'l 


1 7 99('92'l 


15 12 ('66') 




[xlO-2 GeV2] 


12.72(13) 


12.542(63) 


12.522(53) 


12.511(45) 


7 

KPn) 

A{m^,Pn) 


[degrees] 


1.12865(64) 
-17.4(14) 

-0.568(48) 


1.07562(21) 
-13.61(65) 

-0.584(27) 


1.05469(10) 
-12.56(52) 

-0.636(26) 


1.038631(69) 
-11.68(42) 

-0.705(24) 


M^T^^Pn) 




-0.969(81) 


-1.603(75) 


-2.226(91) 


-3.13(11) 



TABLE XVII. Results for n = 1 state in the laboratory system L2 with energy state 
cut-ofF N = 3. Two scattering amplitudes are defined by yl(m,r,Pn) = i^''^^{Pn)/Pn 'En/2 and 
A{mTr,Pn) = (/i"*//7r)^ ■ ^("^TTjPn)) where /^"* is the pseudoscalar decay constant measured on 
lattice and f^, = 93 MeV. 
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24^ X 48 


1.47 


0.1076(13) 


2.583(31) 


0.1357 


0.806(1) 
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395 












0.1367 


0.755(2) 


4000 


10 


390 












0.1374 


0.691(3) 


4000 


10 


380 












0.1382 


0.576(3) 


4000 


5 
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TABLE I. Simulation parameters. Dagger symbol means that we average two measurements 
on the same configuration, one with the temporal origin located at t = and the other located 
at t = T/2. The lattice spacing a is fixed by the p meson mass at the physical pion mass and 
nip = 768.4 MeV. A^xraj is the number of all trajectories, A'^skip is the number of separation between 
two measurements and A^Meas is the number of the configurations used the measurements. 



CM P =( 0, 0, 0) ra = n = 1 n = 2 n = 3 



Pl,n 


( 0, 0, 0) 


( 1, 0, 0) 


( 1, 1, 0) 


(1, 1,1) 


P2,n 


( 0, 0, 0) 


(-1, 0, 0) 


(-1,-1, 0) 


(-1,-1,-1) 


LI P =( 1, 0, 0) 


n = 


n = 1 


n = 2 


n = 3 


Pl,n 


( 1, 0, 0) 


( 1, 1, 0) 


( 2, 0, 0) 


( 1, 1, 1) 


P2,n 


( 0, 0, 0) 


( 0,-1, 0) 


(-1, 0, 0) 


( 0,-1,-1) 


L2 P =( 1, 1, 0) 


n = 


n= 1 


n = 2 


n = 3 


Pl,n 


( 1, 1, 0) 


( 1, 0, 0) 


( 1, 1, 1) 


( 1, 0, 1) 


P2,n 


( 0, 0, 0) 


( 0, 1, 0) 


( 0, 0,-1) 


( 0, 1,-1) 



TABLE II. Momentum assignment for the source operator iln ^ (t) defined by Eq.(15). Here 
Pi^n is the i-th pion momentum of the n state P is the total momentum of two pions system 
P = Pi,n -I- P2,n in units of 27r/L. 
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/3 F [GeV] 


Hn) Cl 


xVd.o.f. 


ao/m^ [l/GeV2] 


1.80 0.135(20) 


1.74(15) 0.70(77) 


0.04 


-1.04(40) 


1.95 0.1035(59) 


1.199(72) 0.765(85) 


1.55 


-1.76(20) 


2.10 0.0909(24) 


0.963(28) 0.702(29) 


4.78 


-2.26(11) 


TABLE III. Results of fitting for the scattering length ao/m^ 


r obtained with the fitting function 


Eq.(33), where we set /x = 1 GeV. The results at the physical pion mass are also tabulated. 


/? ^00 


^10 [l/GeV2] 


A20 [l/GeV^] 


A /U.U.i. 


1.80 -0.033(89) 


-0.78(29) 


0.37(20) 


0.09 


1.95 -0.184(57) 


-0.63(18) 


0.26(12) 


0.53 


2.10 -0.273(35) 


-0.55(12) 


0.169(83) 


5^ 


TABLE IV. Results of fitting for the scattering length ao/m7r obtained with 


d CIlVcl^tJliL iUIlll 


defined by Eq.(34). 








P Aw [l/GeV2] 


A20 [l/GeV^] ^30 [l/GeV6] 


xVd.o.f. 


ao/m^ [l/GeV2] 


1.80 -1.01(29) 


0.79(76) -0.20(46) 


0.02 


-0.99(28) 


1.95 -1.57(14) 


1.58(33) -0.54(18) 


2.1 


-1.54(13) 


2.10 -1.975(82) 


2.18(20) -0.82(11) 


8.2 


-1.932(78) 


TABLE V. Results of fitting for the scattering length ao/m^ obtained with 


a divergent form 


defined by Eq.(35). The results at the physical pion mass are also tabulated. 




P = 1.80 


P = 1.95 


P 


= 2.10 


ml [GeV2] [GeV] 


ml [GeV^j /^"* [GeV] 


ml [GeV2] 


fi^' [GeV] 


1.128(1) 0.2516(11) 


1.287(1) 0.2190(16) 


1.332(2) 


0.1959(19) 


0.814(1) 0.2279(14) 


0.857(1) 0.1908(15) 


0.896(1) 


0.1739(21) 


0.571(1) 0.1993(11) 


0.573(1) 0.1650(14) 


0.605(2) 


0.1540(22) 


0.238(1) 0.1613(15) 


0.291(2) 0.1381(21) 


0.291(1) 


0.1214(19) 


{mfy')^ 0.1287(33) 


{mfy')^ 0.1054(47) 


{mf^y^f 


0.0895(45) 


chiral limit 0.1260(31) 


chiral limit 0.1032(47) 


chiral limit 


0.0869(46) 



TABLE VI. Pseudoscalar decay constant /^"* measured on lattice [16]. Here m^y^ is the 
physical pion mass. 
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/3 


ao/m^ [l/GeV2] 


xVd.o.f. 


1.80 


-2.83(13) 


0.27 


1.95 


-2.543(73) 


0.92 


2.10 


-2.449(40) 


0.54 



TABLE VII. Results of a constant fit for clq = {fi^^/ fn)^ ■ aQ/m^. The results at the physical 
pion mass are also tabulated. 





A 


B [GeV] 


xVd.o.f. 


aomTT 
ChPT 


-0.0558(56) 
-0.0413(28) 

-0.0144(10) 


0.0328(86) 
-0.0119(43) 


0.02 
0.65 


TABLE VIII. Results of the continuum extrapolations for the two scattering lengths 
flo^TT = (/i"*//7r)^ ■ CLoiriTr as a function A + aB. Here A corresponds to the scattering 
the continuum limit. 


aorTT-TT and 
; length in 




1.80 


1.95 


2.10 


^10 [l/GeV2] 
A20 [l/GeV^] 
^30 [l/GeV6] 
^01 [l/GeV2] 
An [l/GeV4] 
A21 [l/GeV6] 
xVd.o.f. 


-1.33(21) 
1.62(53) 

-0.69(31) 

-0.83(44) 
1.4(12) 

-0.46(83) 
0.90 


-1.52(12) 
1.51(29) 

-0.52(16) 

-1.18(47) 
1.9(11) 

-0.65(61) 
0.64 


-1.899(84) 
2.00(20) 

-0.73(11) 

-1.43(40) 
2.48(95) 

-1.04(52) 
1.33 


TABLE IX. Results of a polynomial fit of 
A{mT^,p) = t&n S{p)/p ■ E/2. 


and p'^ for the scattering amplitude defined by 


P 


1.80 


1.95 


2.10 


lio [l/GeV2] 
ioi [l/GeV2] 
in [l/GeV^] 
xVd-o.f. 


-2.84(10) 
-2.78(59) 
3.67(76) 
0.77 


-2.546(59) 
-2.84(44) 
3.36(49) 
0.50 


-2.438(38) 
-2.14(37) 
2.05(43) 
0.75 



TABLE X. Results of a polynomial fit of rrv^ and p"^ for the normalized scattering amplitude 
defined by A{m^,p) = {ff/Uf ■ te.n6{p)/p ■ E/2. 
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p2 [GeV2] 


[GeV] 


S{p) [degrees] 


0.020 


0.40 


-3.50(64) 


0.072 


0.60 


-9.5(30) 


0.140 


0.80 


-16.9(64) 


0.232 


1.00 


-25(10) 



TABLE XI. Scattering phase shift 6{p) in the continuum limit at the physical pion mass. 



27 



FIGURES 



two-pion energy [GeV] 



2.2 



4m 

n 



■ measured 
□ calculated 



2.0 - 



%n (1300) 



1.8 - 



B 



□ 



1.6 - 



B 



1.4 - 



1.2 - 



1.0 



CM 



L1 



L2 



FIG. 1. Center of mass energy of two-pion system for m-jr/nip f» 0.6 at /3 = 2.10 without the 
two-pion interaction. We measure the scattering phase shifts at the energies referred by filled 
symbols. States denoted by open symbols are used only to examine the effect of the cut-off of the 
state for the diagonalization. Solid and dashed lines denote the 7r7r*(1300) state energy and the 
inelastic scattering limit, respectively. 
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FIG. 2. Examples of the pion four-point function Gnm{t) in the center of mass system CM and 
two laboratory systems LI and L2 for mT^/mp ~ 0.6 at /3 = 2.10. Filled and open symbols indicate 
positive and negative value. 
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FIG. 3. Ratios Rn{t) and Dn{t) for n = in the center of mass system CM with energy state 
cut-ofF N = 2. mTr/nip increases from top to bottom, while /3 increases from left to right. 
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FIG. 4. Ratio Rn{t) and Dn{t) for n = 1 in the center of mass system CM with energy state 
cut-off N = 1,2. rriTr/mp increases from top to bottom, while /3 increases from left to right. 
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FIG. 5. Ratio Rn{t) and Dn{t) for n = state in the laboratory system LI with energy state 
cut-ofF N = 3. mTr/nip increases from top to bottom, while /3 increases from left to right. 
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FIG. 6. Ratio Rn{t) and Dn{t) for n = 1 state in the laboratory system LI with energy state 
cut-ofF TV = 1, 2, and 3. rriT^/mp increases from top to bottom, while (3 increases from left to right. 
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FIG. 7. Ratio Rn{t) and Dn{t) for n = state in the laboratory system L2 with energy state 
cut-ofF N = 3. mTr/nip increases from top to bottom, while /3 increases from left to right. 
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FIG. 8. Ratio Rn{t) and Dn{t) for n = 1 state in the laboratory system L2 with energy state 
cut-ofF N = 1, 2, and 3. m-n/mp increases from top to bottom, while /3 increases from left to right. 
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FIG. 9. Ratios ii„(t) and Dn{t) for n = and 1 states in the laboratory system L2 with energy 
state cut-off A'^ = 3 at /3 = 2.10 for rriT^/mp k, 0.6. 
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FIG. 10. Results for ao/rriT^ at each (3 and polynomial fits (dashed lines). 
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FIG. 11. Results for ao/rriT^ = (/^"V/tt)^ ' (^o/'m-K with /^r = 93 MeV at each (5 and constant fits 
(dashed hnes). 
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FIG. 12. Results for aorriT^ and aom.,^ = (/i"*//7r)^'ao/?^^7r at the physical pion mass as functions 
of lattice spacing. 
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FIG. 13. Scattering amplitude A{mT^,p) = taxi.6{p)/p ■ E/2. CM„ refers to the amplitude 
obtained from the n-th state in the center of mass system, and Ll„ and L2„ to those in the 
laboratory systems. Open symbols are excluded from the global fit. m.j^/mp increases from top to 
bottom, while P increases from left to right. 
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FIG. 14. Scattering amplitudes A(m^,p) = tan5(^)/p-£;/2 and ^^(m^rjP) = {f-^^ I f -nf ■ -^[m-K ,V) 
at the physical pion mass and several momenta. The pseudoscalar decay constant /|"* is measured 
on lattice and /jr = 93 MeV. Open symbols are aom^ and a^niT^ taken from Fig. 12. 
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FIG. 15. Results of the scattering phase shifts in the continuum hmit (dashed Hne) and a band 
of error bars. Solid line [28] is estimated with experimental inputs using the Roy equation. Symbols 
represent data of Aachen- Cern-Munich Collaboration [30] and of Losty et al. [31]. 
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